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Abstract 

The usual derivative expansion of gravity duals of charged fluid dynamics is known 
to break down in the zero temperature limit. In this case, the fluid-gravity duality is 
not understood precisely. We explore this problem for a zero temperature charged fluid 
driven by a low frequency, small amplitude and spatially homogeneous external force. 
In the gravity dual, this corresponds to time dependent boundary value of the dilaton. 
We calculate the bulk solution for the dilaton and the leading backreaction to the metric 
and the gauge fields using the modified low frequency expansion of The resulting 
solutions are regular everywhere, establishing fiuid-gravity duality to this order. 
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1 Introduction 



The AdS/CFT correspondence has provided useful insight into the dynamics of strongly cou- 
pled quantum field theories, particularly nonabelian gauge theories. Recently this has led to a 
fluid-gravity correspondence which provides a study of conformal fluid dynamics, an effective 
description of strongly coupled conformal field theory at long wavelengths in local equilib- 
rium. In fact there is a precise mathematical connection between a long distance limit of the 
Einstein gravity and its holographic dual fluid dynamics [11 [21 El IH El E]- Explicitly, it has 
been demonstrated that certain deformations of asymptotically AdS^+i black branes which are 
slowly varying along the boundary directions (but can have 0(1) amphtudes) provide dual 
descriptions of solutions of the equations of fluid dynamics. 

In fluid dynamics, there are local thermodynamic quantities such as local temperature, 
chemical potentials, i?-charges or d-velocity, which are slowly varying along the boundary 
directions compared to effective equilibrium length scale of the fluid, the mean free path Imfp- 
Bulk dual solutions of this inhomogeneous fluids are constructed order by order in derivatives 
respect to boundary coordinate. 

More precisely, the asymptotically AdS space is foliated into a collection of tubes, each 
characterized by a value of the boundary coordinate. Each tube is centered about a radial 
ingoing null geodesic starting from AdS boundary. The width of each tube in the boundary 
direction is smaller than the scale of dual fluid dynamics. The gravity solution is developed 
locally in each tube, which turns out to be black brane with local thermodynamic quantities 
of the boundary fluid. The global geometry is constructed by gluing the tubes, in which the 
local thermodynamic quantities change along the boundary direction. 

One crucial aspect of this construction is ultra-locality. The local thermodynamic quantities 
and metric corrections are expanded around a point on the boundary, which may be chosen to 
be ^^^=0. For some local thermodynamic quantity gj(x^), 

q,ix>') = g,(0) + x^d^q^iO) + ^x^^x^d^dM^) + ■■■ (1) 

In equilibrium these quantities are independent of x^ so that all the higher terms vanish. The 
bulk then corresponds to some static black brane. These quantities appear in the bulk metric 
and other fields, generically denoted by g{r,x^), where r is the AdS radial coordinate. One 
then expands 

g{r,xf') = g^°\r,xf') + c/W(r,x'') + g^'^\r,x'') + (2) 

where g^^\r, x'^) is the bulk field obtained by replacing Qi in the equilibrium solution by qi{x^). 
g^^\r,x^^) is clearly not a solution of the bulk equations of motion. The higher order terms 
g^"'\r,x^) are then determined by requiring that the full g{r,x'^) solves the bulk equations of 
motion. The expansion above then constitutes a derivative expansion of a bulk solution. The 
equation satisfied by 

gin) is 

of the form 

H{g^'\q,m)g^-\r,xn = s^. (3) 

if is a linear differential operator in the second order in the radial variable. s„ is a source term 
from mth order corrections where m < n. The operator does not contain any derivative in the 
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boundary directions, since this would produce terms of higher order. In other words, in this 
derivative expansion, H becomes ultra-local operator. The corrections g^^\r,x'^) are required 
to be regular everywhere outside the black brane horizon. They are also required to fall of 
sufficiently fast at the boundary r = cxd, so that they do not lead to additional sources in the 
boundary theory. 

In [3lll], this program has been carried out for conformal fluids carrying a global R charge. 
However, it was found that in the tubes where the local temperature is zero, i.e. the black 
hole locally looks like an extremal black hole, the solutions are singular at the horizon. These 
singularities seem to be genuine in that those possibly cause singularities in the curvature 
invariants at black brane horizon. 

The failure of the derivative expansion method of deriving the bulk solution in these tubes 
does not agree with expectations from the dual fluid dynamics. In the dual fluid dynamics, 
the effective equilibrium length scale is the mean free path, which is given by /^/p ~ - [3 E]- 
r] is shear viscosity and e is energy density. For the fluids described by the gravity solutions in 
[U El El in El El [ID], ?7 = where s is entropy density. Consequently, Imfp ~ if^- ^ tube 
where the local geometry becomes extremal, the zeroth order gravity solutions in the derivative 
expansion in these papers have finite entropy and energy densities. This implies that there 
should be a reasonable fiuid dynamics even at zero temperature. 

These divergences appear even in the linearized regime, i.e small amplitude and small 
frequency perturbations around an extremal black brane. Consider for example a linearized 
scalar field in the extremal background. It has been shown in |TT] that in this case performing 
a small frequency expansion in the scalar field equation is not straight-forward. Technically, 
this is because the equation contains terms which have powers of the frequency multiplied by 
functions which blow up at the horizon. This may be traced to the fact that the components of 
the background metric have double zeroes or poles at the horizon as opposed to single zeroes 
or poles for a geometry at finite temperature. 

A clue to resolve this problem is obtained by observing that the near horizon geometry of 
the charged black brane becomes AdS2 x R^'^, where d+1 becomes bulk space-time dimension. 
For the near horizon region, it is natural to introduce a new radial coordinate C, ~ ^r^, where 
oj is frequency of the scalar field, r is radial coordinate of the black brane, and Tq is the horizon 
of the black brane. Motivated by the fact that in the near horizon region, the equations 
involve C with no further u dependence, the small frequency expansion is then obtained by 
writing everything in terms of C and then expanding in powers of u. To lowest order in this 
modified small frequency expansion, the equations are then solved in two regions : (i) the inner 
region close to the horizon, which is defined as C ~^ oo with ^ — )■ 0, and (ii) the outer region 
which corresponds to small Q oi r ^ tq. It is essential to use the coordinate C in the inner 
region rather than r. The solution in the full space is then obtained by matching the inner 
and outer region solutions. By construction, a solution which is regular at the horizon exists. 
This technique has been used to obtained various response functions by solving the linearized 
gravity equations and gauge field equations [121 [13]. 

As will be clarified below, this modified low frequency expansion implies that the tubewise 
approximation used in obtaining gravity duals of boundary fiuid fiows breaks down. This is 
because the change of radial coordinates from r to C involves the frequency, so that in position 
space this implies a non-local (on the boundary) redefinition of fields. This reorganizes the 
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low frequency expansion in which differential operators in the equations are not ultra-local any 
longer. The implications of this fact for fluid mechanics are not clear at the moment. 

In this paper we consider a related problem where similar divergences appear and show 
how these get resolved. We consider four dimensional Einstein-Maxwell-dilaton theory with a 
negative cosmological constant. All bulk fields are spatially homogeneous but time dependent. 
The dilaton has a nonzero boundary value and the velocity of the boundary fluid is zero. First 
we consider slowly varying deformations of a charged black hole in this geometry in the presence 
of the dilaton source. These deformations are entirely due to the dilaton source. The setup is 
similar to that of [13] where a time dependent boundary value of the dilaton evolves an initially 
pure AdS geometry into that with a space-like region of large curvature. In the sense of fluid 
dynamics, our set up is a natural extension of |2j to /^-charged fluid with vanishing velocities. 
The authors in [2] solve the Einstein-dilaton system with negative cosmological constant in 
which the boundary value of the dilaton field is slowly varying with arbitrary large amplitude. 
The field theory dual of the gravity system becomes a certain fluid dynamics satisfying Navier 
Stokes equations with dilaton dependent forcing term. We first treat the problem in a naive 
derivative expansion. In this expansion, only the dilaton can have 0(1) changes, as in [T3]. 
The changes of the metric and gauge fields are due to the backreaction of the dilaton, and 
therefore suppressed by powers of the frequency. In [H] (as in IH [2]), the equations which 
determine the corrections to the fields at any given order n are linear and contain fields of 
order {n — 1). As expected, we find that the derivative expansion breaks down and solutions 
are singular at the horizon. 

We then explore if the modified low energy limit can tame these divergences, keeping the 
variation of the dilaton field 0(1). We find that the effect of backreaction cannot be ignored 
in solving for the corrections to the dilaton and other fields even in the lowest order, even 
though the backreaction is small. This is because the nonlocal change of the radial coordinate 
typically implies that radial derivatives are large (in terms of the parameter of the derivative 
expansion) rather than 0(1), thus making the effect of the backreaction large. 

This motivates us to consider the problem for the case where the dilaton has both small 
frequency as well as a small amplitude. We calculate the bulk dilaton and its backreaction 
to the metric and gauge fields to the leading order in amplitude. We find that in this case 
the modified low frequency expansion of [11] and the matching procedure indeed leads to bulk 
solutions which are smooth everywhere. This setup is similar to [HI |12l [13] . In these papers, 
the linearized problem was solved for fields which also depend on the spatial coordinates on 
the boundary, but the backreaction of the fields were not calculated. In our problem there is 
no such spatial dependence; the change of the metric and the gauge fields are entirely due to 
the backreaction of the dilaton. A consistent treatment of the backreaction however require 
us to go to higher orders in the frequency(in some case up to the fourth order). At the same 
time, we need to go to the second order in the amplitude of the boundary dilaton. The fact 
that this scheme works to lowest non-trivial order in the presence of back-reaction indicates 
that there is a systematic double expansion in the frequency and the amplitude which leads to 
non-singular solutions. 
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2 Charged Black Brane with Dilaton Field 



In this section, we define our model to address the problem. We consider 4-dimensional 
Maxwell- Gravity theory with time dependent dilaton. The solutions are constructed order 
by order in small frequency in the perturbation theory. It turns out that leading corrections in 
the perturbation theory are divergent as they approach the black brane horizon in the extremal 
limit. 

2.1 Derivative Expansion 

A consistently truncated theory from M-theory with S"^ compactification [TB] motivates us to 
consider Einstein-Maxwell-dilaton theory with negative cosmo logical constant, 

S=^J d'x^g (ii? - \fmnF'^^ + ^ - ^^M^a^V) , (4) 

where /€4 is the gravitational constant. Indices M,N.. run from to 3, and Fmn is field 
strength from U{1) gauge field Am- We choose units with L = 1. The equations of motion are 

Wmn = Rain + ^Omn — '^FmpF^ + -QmnFpqF^'^ (5) 

- ^dM4>dN4> = 0, 

= VmF^'^ = 0, (6) 

X = VV = 0, (7) 

where V denotes a covariant derivative with bulk metric Qmn- A charged black brane solution 
of these equations of motion in Eddington-Finkelstein coordinates with constant dilaton is 

- Uo{r)dv^ + r^dx'dx\ (8) 

- I) dv, (9) 

(10) 

where UQ{r) = + ^ — y and tq is the outer horizon of the black brane, the largest root of 
Uo{ro) = 0. V is the ingoing null coordinate which is time coordinate in the AdS4^ boundary, p 
and e are charge density and energy density of the black brane respectively. 

Let us consider time dependent dilaton which is slowly varying compared to tq, (j) = 0(o)(f). 
More precisely, the dialton has a form of 

0(o)(t^) = /(-), (11) 
where e is dimensionless small parameter. The function / satisfies 

/'(-)- 0(1), (12) 



ds^ 


= 2dvdr 


A 












= const, 
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where prime denotes derivative with respect to its argument. The derivative of the dilaton 
with time is suppressed by e. 



V(-)~-~OW. (13) 



dv ro ro ro 

(t>(o){v) is obviously not a solution of the equations of motions. To solve the dilaton equation 
perturbatively, we add correction terms. The dilaton is expanded as 

(f){r,v) = (j)(^o){v) + (f){i){r,v) + (f)(^2){r,v)..., (14) 

where ^(o) is zeroth order in e, and 0(i) is 1st order in e and so on. The dilaton solution can be 
calculated perturbatively order by order in e, which becomes the expansion parameter of the 
perturbation theory. We promote the energy density and the charge density to be functions of 
time as 

p{v) = po + rlC{v), (15) 
e{v) = eo + r^,E{v), (16) 

where po and eo are constants. For further convenience, we set C(— oo) = E{—oo) = as initial 
conditions. We expand C{v) and E{v) as 

C{v) = C(°)(^) + C^'n^) + C^'H^)- (17) 
E{v) = E^''\v) + E^'\v) + E^^\v)... (18) 

The dilaton equation up to first order in e becomes 

= dr {r^U{r,v)dr<Pii)ir,v)) + 2r9,0(o)(t;), (19) 

where U{r, v) = + (Po+''of("))' _ ^i^o+jMuR, As we will show below (See Eq([35D and EqdSSD), 
C{v) and E{v) are higher order in e. The first order correction to the dilaton is given by 

<l>ii)ir,v) = r "'^li'!\~"\ d.<Pio){v))+A2{v), (20) 



r'^Uo{r) 

2 

where Uq^t) = r^ + ^ — ^ and Ai(^;) and A2{v) are integration constants which are to be 
determined by boundary conditions that we demand. The regularity condition at the black 
brane horizon requires Ai{v) = 1. Moreover, we want a specific boundary condition that as 
r — )■ cxD, (f){r,v) = 0(o)(f). A2(t') is determined by this boundary condition. 

We pause here to demonstrate the relationship of the tube-wise solution with derivative 
expansion. Consider the congruence of null geodesies emanating from AdS^^ boundary and 
tubes which are centered along the null geodesies. The set up is spatially homogeneous, so we 
classify the tubes by v. Without loss of generality, we set = for every individual tube. We 
expand 0(o) {v) in the neighborhoods of = as 

0(o)(^^) = 0(o)(O) +£t;a,(/)(o)(0) + l£Va20(o)(O) + ... (21) 
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The charge density and energy density can be expanded as 

p{v) = p(0) + evd,p{0) + ^e'v'dlp{0)..., (22) 
e{v) = e{0) + evdye{0) + ^e\^dle{0)... (23) 

We add correction terms to the dialton field as 

^(r) = 0(0) + £0(i)(r) + £'0(2)(r)... (24) 

We omit w-dependence in the correction terms because the differential operator acting on 
these becomes ultra- local as argued below Eq([3]). Plugging these into the dialton equation and 
evaluating it up to first order in e at f = 0, we obtain 

= dr {r^U{r)dr(Pii) (r)) + 2r9,0(o) (0) , (25) 
where U{r) = + — The solution of this equation becomes 

0(1) W = J \%[^^ (g.0(o)(O))+A2(O), (26) 

where Ai(0) = 1 for the regularity of the solution 0(i)(r) at the horizon and A2(0) is determined 
by the same boundary condition of our solution at AdS^ boundary. To get global solution, we 
should patch every local solution. Even if the way of getting solution is different, Eqf l20|) and 
Eq( l26i) have the same form at least in the first order in small frequency expansion. This is also 
true for its back reactions. 

Back reactions are obtained perturbatively order by order in e with gauge field and metric 
being expanded as 

_ (0) I (1) I (2) /o'7^ 
9mn — 9mn + 9mn + 9mn--- l^'J 

Am = + 4^+4?... (28) 
For leading order correction, we try following form of metric and gauge field solutions: 

ds'^ = 2dvdr - U{r, v)dv'^ + r'^dx'dx' + ^^-^ ^ dv^ - 2h{r, v)drdv, (29) 

A = p{v)( — -— )dv + a{r,v)dv, 

ro(f) r 

where k{r,v), h{r,v), and a{r,v) are leading order corrections in the perturbation theory. 
Details of equations of motions and calculations of the solutions are in Appendix El We briefiy 
list the leading back reactions. The leading corrections to gauge field and metric are 

9i^N = <^ = 0, (30) 
Hr,v) ^ = -\id^Mi^)r (^^^^ dr' + h{v), (31) 



/r ( 2 i2\2 

+ 02(1;). 

ai(t'), (32 ('*^) and ki{v) are integration constants. They are determined by specific bound- 
ary conditions that we demand. At the black brane horizon, these solutions are already regular 
by choosing Ai(f) = 1. For the boundary condition at r = 00 we demand that each leading 
correction of the perturbation theory behaves as 

h{r,v) ~ 0(r°), (34) 
k{r,v) ~ O(r^), 
a(r, t>) ~ 0(r^^). 

The motivation behind these boundary conditions is that there are no non-normalizable modes 
which deform the boundary metric, chemical potential or charge density [H [21 |3l Uj. 

There are the constraint equations which are certain combinations of the equations of back 
reactions. They are in fact the equations of dual fluid dynamics [2]. 

C(o)(t;) = C^^\v) = E^''\v) = 0, (35) 
E^'\v) = ^{d.<l>^o){v)r. (36) 

The other components of gauge field and metric are trivial. 

2.2 Divergences of the Leading Order Corrections in the Extremal 
Limit 

The regularity of the solution that we impose for each perturbative correction in the previous 
section breaks down in the background of extremal black brane. In Appendix [B], we derive 
the near horizon behavior of leading corrections of the dilaton and its back reactions in the 
extremal limit by expansion in u — 1, where u is a rescaled radial coordinate, u = ^. To be 
more general, we keep Ai(f) to be arbitrary in Eq fll39p . As shown above, Ai{v) = 1 ensures 
regularity at the horizon for non-extremal black brane. However in the extremal limit, all the 
corrections, Eq (l20|) .Eq (l3T|) .Eq (l32!) and Eq (l33|) have singularities at the horizon. For example, 

ct>{u, v) = 0(0) - ?^^^^ln{u - 1) + 0(1), (37) 
org 

as M — )■ 1. The near horizon expansion of the back reactions are also leading to divergences 
in physical quantities like curvature invariants and field strengths. The m — )■ 1 behavior of the 
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leading correction to the gauge field is given by 

a(u, v) = _^ (9-ho)i^)y /g^^/ _ ^^ _ _ _ 1) + o(u - 1)2) + 0(1). (38) 

108 ro 

The first two terms can cause singularity in the field strength, Fry The divergences of metric 
corrections are 

1 (a,0(o)(t;))' f 1,2 

k{%v) = -\rlid.<l^^o){v)r(^liu-irin{u-l) + 0{l)y (40) 

In particular the term multiplying (u — iyin{u — 1) in k{u, v) can possibly cause singularities 
in curvature invariants. 



2.3 Divergence Resolution (The Main Result) 

In this subsection we briefiy discuss the main result of this paper without much technical 
details. To deal with divergences discussed in Sec J2.2[ we follow Ref. [TT] and divide the radial 
coordinate into two regions. 

Inner Region : w — 1 = — for (5 < < oo, (41) 

Outer Region : — < n — 1, 

with a certain scaling limit, 

z/ ^ 0, C = finite, 5^0, and ^ ^ 0. (42) 



where u is frequency of the fields in the perturbation theory. Note that switching the radial 
variable m to ^ is a non-local transformation. For this transformation, we need to evaluate 
our equations in the frequency space. We use ^ as a radial coordinate for the inner region 
and u as that for the outer region. We also define overlapping region (or matching region) as a 
region that ~ 5. The black brane horizon is located at = oo. As — )■ 0, we approach the 
overlapping region by the scaling limit f H2|) . The solutions listed in Sec J2.1l are the outer region 
solutions. The expansions of the dilaton and its back reactions as u — 1 in Sec 12. 2 1 are therefore 
the fields in the overlapping region. With the scaling limit we define a new perturbation theory 
in small frequency in the inner region. If the inner region solutions are 

• Regular at the black brane horizon, 

• Smoothly connected to the outer region solutions in the overlapping region, 
the solutions are regular everywhere. 
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The dilaton equation in the extremal background is 

1 / 1 \ 

VV(«, v) = ^{ duiu^Uiu, v)ducl)iu, v)) + -d^iu^d,(l)iu, v)) + —d^d,<j)iu, v))=0, (43) 

where again we define dimensionless radial coordinate, u = The metric factor U{u,v) = 
^"^2^-* [u^ + 2u + 3) + ^^'"*'^2^'" ~ ^^T^- More explicitly, the equation has a form of 

1 

= du({u-l)\u^ + 2u + 3)du(j){u,v)) + -du{u^d,(f){u,v)) + —dud,(j){u,v) (44) 

ro To 

- 2E{v)du{udu4>{u,v)) . 

In this equation, we see that the term multiplying the energy density E{v) is of 0{e'^) whereas 
the terms proportional to 5^0(u, ti) is of 0{e), using the constraint equation fl5^ (which shows 
E{v) ~ 0(e)) and the dilaton solution f l2U]) . Then, it may appear that the last term in Eq flU]) 
can be ignored for obtaining the first order solution in e. However, this is no longer true when 
we switch the radial variable u to ^. The momentum u which appears in Eq( HT|) is effectively 
proportional to e. This is because the dilaton field is localized around ~ e in the momentum 
space(See the discussion in the beginning of SecJl]). In ^ coordinate, each it-derivative in Eq (Hl]l 
produces extra factor of -. Therefore, the first term in Eq (H^ becomes the higher order in e 
than the term multiplying the energy density in C, coordinate. This means that this later term 
cannot be ignored any longer. We have not yet been able to solve this type of equation. In 
the following we will choose a regime where the amplitude of the dilaton field is small. This 
will allow us to ignore the last term in Eq( H^ . but retain the essential feature of the problem. 

The inner region solutions that we solve are completely agree with above two conditions 
for the entire solutions to be regular everywhere. To be precise, we briefiy list our inner region 
solutions. The inner region solution of the dilaton has a form of 

0M.(O = 0t°^ + ^(A«-^e^«i?i(^o)+^^4^^ (45) 

J oo ^ J 

where 




, , liAf^ ^id)f^ lid)]^^ ( 4zA , if , 

ro 3ro4 3ro^ V 3ro/ 3ro 
Ek{x) = J-|°° ^jr-dt is The Integral Exponential Function, where k is an integer. The numerical 
coefficients, aI}\ A^if^ and so on, are determined by matching with the outer region dilaton 
solution ( 137|) in the overlapping region, more explicitly with Eqf ll39p . The inner region solutions 
are developed in the frequency space. Therefore, for matching we need to take the outer 
region solution to momentum space by Fourier transformation defined in Eq (!59i) . The precise 
expression of the outer region solution of the dilaton in frequency space is in Sec J3.2[ We have 
a precise form of A^u'^ as 
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where 7 = 0.57721.. is Euler's constant. For the regularity at the horizon and the matching 
with the outer region dilaton solution, only ingoing waves at the black brane horizon are 
allowed in Eq fHSl) . As a homogeneous solution of the dilaton equation we could add outgoing 

waves as M^^e^-o to EqdlS]), where S^"^ is a constant of order n in e. The matching condition 
forces B^^ = 0. For n > 1, it turns out that bI^^ spoils the regularity of the dilaton solution 
in n + 1th order in e. Thus, we naturally impose ingoing boundary condition at the horizon for 
the smooth dilaton field. The other coefficients are obtained by matching with higher order 
solutions in the outer region. Near horizon, the inner solutions behave as 

0(.n).(O - + + (48) 

+ ^2 

which is manifestly regular as — >■ 00. 

We evaluate back reactions from the dilaton up to the second order in small frequency as 

h^nAO = C^lO + '^^^^^^ + ^i.lO-, (49) 

J 00 y s 

hnUO = + Ul^' + ^)+ - 2 £ pZ^) (51) 



o A 3 (a /"^ 1,(2) / A , Q dy <i) ^ ^fn),.(0 o ^Mz^(^) 
- 2roZ^ [6/ ^/i(4.(?/) + 3y^^/i;,^),(2/)-3^^ 2 ^3 



00 



+ 12ry [ ^ r^h\'\ (.)..., 
^00 y J 00 ^ 



where 



uj{u - uj) y'^ ^_i_y ^ ^ iuy - uj)y. 



u{u - u)y'^ iuy - u})y^ 



Q 2 2 e^oyE,{^)E,[ 
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, ,(0) u'^iv — uY [^0,1 ^ ^ 

duj<pf_^ ^ , ' / y^dy ^e^-^o^- 



3ro 3roz/ / 

are integration constants which are determined by match- 
ing with the outer region solutions. Egf liQ]) . Eqf lSOj) and Eqf lST]) are smoothly connected to 
frequency space expressions of the outer region solutions of Eq fj39|) . Eq fl38l) and Eql HOj) respec- 
tively. Details are discussed in Sec J4.2[ We determine some of the integration constants in the 
inner region solutions with this as 



Hi'' = ^/ du^^^^<P^^'<P^^l{8 + ^-13V2tan-\V^^ (54) 



864ro^ """" 1/2 
+ 16ln{u)) , 



A^? = -^J^ du^^^^^ct>fct>T-J-l^ + 2V2tan-\V-^^ (55) 
- %ln{y)) , 

Kf' = \rl jyu^-^^<^f<^^l[V2t^^^^^ (56) 

A^^\ Ri'^ and Ri'^ are determined by higher orders in the outer region solutions, is the 
Fourier transform of hiv). h!)^\ hI)'\ A^u\ Af\ kP and X^^^ are finite numerical constants, 
which cannot be obtained analytically(In principle we can. See Sec J4. II for the discussion about 
these constants). The near horizon behaviors of the back reactions are given by 

%„).(0 ~ uHl'^ + u'Hi'' + u'Hi'\.. + 0{^), (57) 

ai.n)AO - u'A^^' + iy'A'^^ + j.'A(^' + u'A(^\.. + 0{-^), 

and %„).(0 ~ u'Rl'' + u'Rl'' + u'ki'' + u'Ri'^ + u'Kl'\.. + 0{-^^^ 

Thus the solutions are regular everywhere. 

Finally, the constraint equations in the inner region turns out to be the same with those in 
the outer region (!35|) . which have forms of 

C^f = C«=i?f = 0, (58) 
1 

E^^ = / duuj(u — u)(j)''^'(l)^^l^ and so on..., 

4^ro y_oo 



in the frequency space. 
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3 Divergence Resolution of Dilaton Field 



In this section, we solve the dialton equation in the inner region. The scahng hmit( l42l) cannot 
be apphed to Eq (H^ because switching radial variable m to ^ is non-local transformation. To 
deal with this, we need to rewrite the equation in the frequency space by Fourier transformation 
as 

(j){u,v) = / e^'^>^(M)rfa;, (59) 



oo 



where (j)uj{u) is a localized and normalizable function in frequency space. For example, we can 
choose (j)u]{u) to be 

(PMr-.e-^e-^f{j)gM, (60) 

where guj{u) is a function carrying radial variable u and f{j) is arbitrary 0(1) function of in 
frequency space. Then, (f}^{u) is suppressed as u approach either zero or oo by the exponential 
factors in it. Consequently, this function is extremely localized around u = ±£ and normaliz- 
able. Using the argument in Appendix [Dl one can show that Fourier transformation of Eq fl60|) 
becomes a form as Eq ffTTj) in the frequency space. This shows that the properties of Eq fl60|) 
is consistent with those of ^(o)(i') introduced in Sec j2.1[ The dilaton equation in momentum 
space can be read off by acting an integral operator ^ e~^^'"dv on Eq fj44|) . The equation 
in the momentum space has a form of 

2 

= du{iu-lf{u^ + 2u + 3)duMu)) + -duiu^Mu)) + duM'^) (61) 



oo 



2 / du {udu(pu;{u)) Ey^^du, 



oo 



where E^,-^ is Fourier transfor of E{v) as defined in Eqf llSSp . In the inner and outer region, 
we expand the dilaton field as 

Inner Region: 0(,„).(O = ^^.(O + + + (62) 

Outer Region: (j)^{u) = (p^^^ + iy(p^J\u) + u'^cp^^^u) + (63) 

respectively. To avoid confusion, we do not tag outer region solutions with ''{out)" . 



3.1 Inner Solution 

As discussed in Sec J2.3l we solve linear dilaton equation by ignoring the last term in Eqf l6ip . 
Switching radial variable u to ^, the dilaton equation in the inner region becomes 

eo^ ((6 + y + ^)5c</'M^(0 j - ^(1 + -^)%<Pi^n)uiO + ^(^ + ^)hn)AO = 0. (64) 

The zeroth order equation in ly is 

Ge^<I^M) - ^-^9,41^.(0 = 0. (65) 
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This equation gives two linearly independent solutions, 

<^t> = ^^"^ + (66) 

where A^^^ is purely incoming wave and the term multiplying B^^ is purely outgoing wave at 
the extremal black brane horizon: ^ = 00. The subscript 'V denotes that the integration 
constants depend on frequency. These two independent solutions are regular at the horizon up 
to the zeroth order in u. First order equation in u is 

t2 



- ^5.0S2).(O + ^ {at + Bl%^^) = 0. (67) 
The solution of this equation is 



^(1) 



(6^ 



where Al and BI are integration constants, which are corresponding to incoming and out- 
going waves at the horizon respectively. Using 



X px' —lax" 

e*""'dx' / —dx" = e"'''EAiax), (69) 



Eq (!68i) becomes 



•/i(o) ■ -0(0) 

^.iO = AT - ATbI%^ - ^e^^i?,(f - (TO) 

^ ' 3ro 3ro Stq Sr, 







The first and the third terms are incoming waves at the horizon, whereas the others are outgoing 
waves. Let us discuss the asymptotic form of the solution. For large y, Ei{y) is expanded as 



EM = —Y.-^^, (71) 
y ^ — ^ y^ 

^ n=0 ^ 

where y is pure imaginary number, so term is bounded. It is manifest that Ei{y) is 
regular as y approaches infinity. In the case that y goes to zero, the function Ei{y) becomes 
divergent. Let us argue what the leading divergence is. To see the leading divergence, we 
calculate following object: 

dEijy) 

lim ^^^^ = lim '^^ = -1 (72) 

where we have used Hospital's theorem for the first equality. Then, the leading divergent term is 
logarithmic. Consequently, this solution is regular at the horizon and divergent logarithmically 
near the matching region. 
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3.2 Outer Solution 

In Appendix [HI we evaluate the dilaton in outer region as m — )■ 1 in the extremal limit. The 
solution in momentum space is given by 

Mu) - + = 4°^ + r ^^^,du + A„ (73) 

To J U^Uq[u) 

where Ai and A2 come from the Fourier transforms of the integration constants Ai(t;) and 
A2(f). They depend on v. Near matching region expansion of Eq fl73l) is 

Uu) = + '""^^ , (1 - Ai) - '-^Hu - 1) (2 + AO (74) 

■ ,(0) / N 

+ (Ai - 7)V2tan-\V2) + 2(Ai + 2)ln{6) + ^(7 - Ai) 

36ro V V2 / 

+ ^(n-l) + 0(n-l)^. 
18ro 

3.3 Matching 

For the inner solution to match the outer one, we need to switch the radial coordinate ^ to u 
near matching region. The matching region is defined as a region with ~ 6. The scaling 
limit fH^ shows that can become a small expansion parameter near matching region. The 
outer region solution is perturbative solution order by order in u. It is justified that one can 
do series expansion of (p(in)u in for matching the outer solution. The inner region solution up 
to the leading order correction in Eq (!70|) . in the radial variable u is given by 

(pi^n)u{u) = Al^^ + Bi°^e^^^ + u(^A[}^ -A^^^Bl'^e^^^ (75) 



.e3ro(u-i)E^{ -) Ei{: 



3ro 3ro(M-l) 3ro 3ro(u - 1) 

Using asymptotic expansion of Ei{y) for small y as 

-i)v 



E,{y) = -7 - ln{y) - (76) 

' ' Tint 



nn\ 

n=l 

we expand the inner region solution in terms of ^j^. The expansion has a form of 

hn)Au) = AW + + f 'f^ - ^{Af + B^^^)ln{u - 1) + (77) 

y3ro(M - 1) 3ro 

15 



We compare the asymptotes of 0(^)1/ (^i) with Egf TM]) to determine each coefficient in it. At 
the zeroth order in 

Af + B^^^ = ct>f. (78) 

At the first order, 

+ = (79) 

B^^^ = (80) 

= ^(^(Ai-7)v^tan-i(v^) + 2(Ai + 2)M6) + ^(7-A0) (81) 

- ^(AW + S(°))-^/nf^V^^4-— 
3ro 3ro \3ro J 3ro V 3ro 

EqdZH]), Eq([79D, Eq([80]) and Eq([8l]) provide 

Ai = l, A(°) = 0(°), S(°) = o, (82) 

and 

4. - ^ f - M.) - ^) . 1^ . f (^) . (83) 

3.4 More on Dilaton Solution in Inner Region 

As is clear from Eq fl77l) . the lowest order solution in the inner region, expressed in terms of 
coordinate u contains all power of ly. This is true for all higher order corrections to the inner 
region solution as well. More precisely, expressed in terms of u will have terms of 

0{iy"^) with m ^ 1. These terms are crucial in ensuring a smooth matching with the outer 
region solution. For example, in Eqf l74l) there is a term ~ u(f)^u\u — 1)- but such a term is 

not present in i^4>\ln)ui'^) ■ show that such a term is actually present in z^^0(^^)j,(m) 

with precisely the correct coefficient. To see this, let us evaluate the second order correction 
of dialton field. The second order equation is 



6^'5|C).(0 - ^-^d,ct>%M) + 9i2)A0 = 0, (84) 



where g(2)u{0 is 



9iM0 = 4e0[2;(O - 4(1 + ^)0S.(O + + (85) 



The prime indicates derivative respect to ^. The solution of this equation is 
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(2) 

where we set Bl ^ = because it cause logarithmic divergence at the horizon in the third order 

(injv 



('21 

in v. Near horizon, (puL^{C) is expanded as 



i7) 



where a° = A^u \ al = —Ai}\ al = — Sir^Aly^ and so on. Regularity condition of 

at the horizon forces bI^'' = 0. Only incoming wave is allowed in ^oo- Near the 

matching region, we switch the radial variable ^ to u for matching with the outer solution. 
Defining a new integral variable y a.s C,' = ^pj, Eq fl86|) becomes 



For matching, we expand 4>(j'\Au) in ly as 



d^) _ 

" j=0 



OO j 

+ E (^7^1)7 (/^]? + - 1) + (^^(^ - 1))') , (89) 



u) 



where /3jJ ,/3jJ and are 0(1) constants, some of which are given by — 

/^o ~ ~§i^ ^'^ '^^^ ^^^^ term in Eq (!89l) is proportional to u, which matches the last 
term in Eq fTT^ precisely. We expect that similar mechanism ensures matching of the higher 
order terms. 



4 Divergence Resolution of Back Reacted Metric and 
Gauge Field 

In this section, we extend our discussion into back reactions. In AppendixjCl we obtain the 
equations of the back reactions without ignoring time derivatives. These equations are the 
starting point for our discussion. 



4.1 Inner Solution 

We begin with Eq fll54p . To solve this equation in the inner region, we need to substitute the 
inner region solution of the dilaton field into it. To do this, we take the inner region dilaton 
solution back to the outer region: (f)(in)uj{0 (t}{in)uj{^^) and plug it into Eq fll54p . With this, 
Eqf ll54p becomes 

duK{u) = -J / dudu(t)(in)y~u{- du(t) {in)u{ ^ -. ) • (90) 
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The w-derivative acting on the dilaton can be switched to derivative with respect to its argument 
as 

u — \ [u — 1)'^ ^ ' u — 1 

After this, we replace the radial coordinate u with ^. Then, Eq( ll54p becomes 

hU.{i) = [J ^ duuiu - J , (92) 

where the prime on the dilaton denotes derivative with respect to its argument, whereas the 
prime on /i(j„)jy(^) does derivative with respect to ^. Plugging the dilaton expansion ( l62l) into 
Eq (l92i) . it becomes 

hUAO = £,(1 + f ) y_ d.u\u - uf (0SJ^O0S2;-J^O (93) 

+ -C.(^oC_(^o + (^ - c.)C.(^?)C-.(^o 
+ ...). 

h{in)u{0 is also localized function around z/ ~ ±£ in frequency space as the dilaton field (See 
the discussion below Eq fl60l) in Secl3]). Then, we expand /i(m)i/(0 as 

h^n)AO = H. + /i£).(0 + ^h%^M) + (94) 

where ^(]^),^(0 is in the first order in small frequency and vh^^in)u^^) second order and so 

on. We expand the other corrections of the back reactions in the same way. Hy is an integration 
constant which can be expanded as = uHu^^ + u'^Hu^K.., where hI) and hI)^\.. are 0(1) 
constants. We note that counting power of e of the solutions to show that each solution is in the 
correct order in small frequency expansion is not manifest in the frequency space. In Appendix 
[D| we discuss details about this power counting by scaling all the frequencies appeared in the 
solutions with e. The solutions up to the second order in small frequency are 

,(1) fC^ _ ,,^7(1)^ f^^^ r ^,.„.,2f, ^2^(1)' 



As will be shown below, the terms containing ^'-integrations appearing in Eq (p5!) and Eq (p6|) 
vanish as — )■ oo. However, the near matching region expansions of those terms present 
additive constant terms. H^"^ and are 0(1) numerical constants which are designed so 
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that they precisely cancel those additive constant terms. With such a choice of h!)^^ and h!)^^ , 
the only constant term in the overlapping region expansion of h{in)v{0 niatch that in the 
outer region solution hy{u) become H„ (See Eq (p^ and Eq fll24p ). There will be numerical 
constants as A^y \ A^u \ Rf^ and R'^f^ appearing in a(^in)u{C) and k(^in)u{C) to be determined by 
the same manner. They provide the near matching region expansions of a{^in)v{,C) and k(^in)u{C) 
to be exactly given by Eq fll26l) and Eq fll27l) . We cannot determine these numerical constants 
analytically but in principle one can obtain the precise values. 

Let us discuss regularity of the solutions. H,^ is regular everywhere. To see the behaviors 
°f ^(2)1.(0 and we expand and in Eq([95]) and Eq([96]) in the hmit of 

large value of their argument x using that Eq(|70l). Eq( |7T]) . Eq llHTl) and B^^ = B^^ = 0. They 
are given by 

J^>' (,) - _ig^ftlM (g,) 



3^0^ n=l 



Ena , , , 



n=l 



Substitution of Eq flWl) to Eq fl^?^ provides near horizon expansion as 



= -A" + E tS^. (99) 

n,m=l 

where 

Amn = l!!^±izlL^,m,n-^ f ^\ ^ <P^^^ cP^^W^ (u - uf^du^ (100) 

4 m + n-l 

where m and n are integers. We note that the cu-integration in Eqf llOOp seems to have poles 
at oj = and u = u for m + n > 2 and could lead to an infinite integrand. However, 

2 

we set ~ as a; —t- as discussed in the beginning of Secj3] (See Eqf l60|) ). This 

ensures that the integration is finite. Then, near horizon behavior of ^(i^)j,(0 is given by 

We obtain near horizon behavior of h^^^-j^iC) by plugging Eq fl97|) and Eq fl98l) into Eq fl96l) . 
which is given by 



n,m=l ^ 



TTi -\- n 1 AjYifi ^ B^^ii^ I (101^ 



m+n tm+n—1 



where 



lm!(n)(-ir+\_„_2 / I V-"" r fAO) „„,,2-n,„ .M-m 



+ 0W<_,a;i-'"(/v-a;)2-")du;, 
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The a" in Eqf fTU^ are proportional to ^L"^ (See Eq([H2D, Eq([HS]), Eq([H7]) and discussion 

(2) f 

{in)u ^ 



below it). This ensures that the integration in Eq fll02p is also finite. Consequently, h(-\{^) 



uHI^^ + 0(|) near horizon. 

Secondly we solve Eqf ll59p which provides solutions of gauge field corrections in the inner 
region. Switching radial variable u to C,, Eqf ll59p becomes 

= ^^T^ + V^%n).(0) • (103) 

We expand the charge density as Ci, = Cl^^ + lyCi^^ + u'^ci^K... a(m)j,(^) can be expanded as 
a(m)i/(0 = Au + + ^'^ul)uiO + ^'^'^u!i)uiO---^ where again is an integration constant 



i exp; 

are given by 



which is also expanded as A,y = vAi^^ + u'^A^u^ ... The solutions up to the third order expansion 



= -roi^^, (104) 



<.g,„K) = -4^' + V3''o£fC«')-4^/fi" + '-o-(^-^j. (105) 



r^(0) ^(1) ^(2)~ 

1-^1/ Ly;/ Oi/ 

- rol^l^^ 7^ + 



The reason why we need to obtain a(j„)z/ (0 to the third order in the small frequency is that 
when we plug the nth order solution of h{^in)u{C) into Eq fllOSp . we get n + 1th order solution of 
a(in)v{i)- For the same reason, we need to get the inner solution of k{in)u{C) up to the fourth 
order in small frequency. 

Let us explore near horizon limit of the solutions. CL~^in)y{C) is manifestly regular at the 

horizon. Using Eq( l99|) and Eq fllOip . we evaluate near horizon expansions of and a(j^);^ 

which have forms of 



o?\ = V. 

(m)v 



m,n=l \ / 



V 1 



- + V3,.„ ± - . V3,.„.|l (108) 



Jji2) /^(O) ^(1) ^(2)X /^(l) ^(2) 



VSroZ/-^ rou ( ^r:^ — + ) + \f?>rQV 
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Then, a(^in)u{0 is finite at the horizon. 

Combining Eqf llSSp and Eqf ll60p . we get 

= -6{u^ - l)K{u) - u{u^ -4u + 3)h'{u) + \ {uk'(u) - K{u)) + 2^3^, (109) 
which gives solutions of k(^in)u{0- Changing radial coordinate u into ^, this equation becomes 



lOu z/3\ /4z/ 6z/2 4z/3 , , , 

= ^'{^ + — + ^ + ^,)hUuiO-Gi^[^j + ^ + ^ + -,)h^n)AO (110) 



r5 



We expand as = + ,(0 + + where is 

a homogeneous solution of Eq flllOp which satisfies 

uKio + ie + ^oKio = 0. (Ill) 

The solution of Eq flllip is given by 

MO = r^Ki^' + (kI^^ + T") + (^^'^ + T") ^^^^^ 

where kI^\ k]^^ and K^^^.. are arbitrary 0(1) constants. We solve Eq flllOp up to fourth 
order in small frequency which are given by 



^(2)1/(0 = 




^(ii)i/(0 = 


2v^roS^^^ 




J oo 



(113) 
(114) 

4 f^KfhW (c'\^^uW ^^0^^,,^'^ 



^ ( hf\ iO - -h'^\ io] + I24u^ (115) 



r<(2) 

2v^roV^, 





/•«' de ( 


J oo C 


L \ 



(116) 



r, 4 /12,(2) /^./N I 6 (1) /^/N _ o;,(2)' _ ^7,(1)' { ti\ 



o-(2) o-(l) ^(3) ^(3) 

+ 12roV^ + 8roV^ - 2v^roV^ + '''^ 
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Kv{C), and are manifestly regular at the horizon. We list behaviors of k^^-l^^iCj 

and k'fK iC) near horizon: 

m,n=l 

r7(2) , rr(2) fr(l) , fr(l) ^(3) ^(3) 

= .K^^) + 12roV ^- +8ro^. ^- + - 2v/3roVi^ + (118) 



4 >^ 6(m + n + 3) .B^^ 4 ^ 



/lOm + lOn + 2 



m,n=l 



+ n + 1 ^rn+n+l U fm+n+2 V ^ + ^ + 2 

m,n=l ^ 



^ 6(m + n + 3) _^ 24(m + n - 1) 



(m + n + 1) (m + n + 2) (m + n) (m + n + 2) 

As ^ — oo, k'^(2i)u^0 ~ i^-^!^^'* + and fc(jti)j,(0 ~ ^^-^i-^'' + 0{j). Therefore, the inner region 
solutions are regular solutions. 

Finally, we solve Eq( ll49p to obtain the constraint equations in the inner region, which are 
given by 

= S(°)-y3C(°), (119) 
-2v^C« + 2^^) + ^j+^J dcoiu - uj)u^f_J^^\ (120) 

-2V3Ci'^ + 2^2) + ^ - / du;{u - uj)uj'<P^:Ia(^\ (121) 



4.2 Matching 

In this subsection, we show that the inner region solutions solved in the previous subsection 
are smoothly connected to the outer region solutions in Sec J2.1[ Let us start with /i(m)i/(0- 
Near matching region, we expand Eq (p5|) with small ^ using Eq (170|) and Eq (!76l) . Integrating 
by ^, Eq (p5|) becomes 

1 / ' \ ^ /*oo 1 

^(H.(^) = 4(^j J rf^^(i^-^)0L°V?i^+ subleading terms..., (122) 

where the "subleading terms" denote terms which are higher order in u when h^(^^-^^{^) is 
expressed in terms of u (We obtain leading corrections only in the outer region). The same 
procedure is applied to Eq(l96l). Near the overlapping region, the expansion of hj^in^^iu) becomes 

''^M-^^) = -I (^) / duu{u-u)(l)l°^(P';^ljn + subleading terms... (123) 
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Combining these, we get near matching region expansion of hi^in)y{u) which is a form of 

(124) 

As M — )• 1, the outer region solution of h^, in momentum space is expanded as 



/^.M = i^^/j I 8 + ^^(125) 

- 13V2tan"^(V2) - 8ln{6) 

(See Appendix [B|. Consequently, Eq 0124p completely matches Eq (]125p requesting that hI^^ 
0, uml'^ = ggi^ J^^ duuj{iy - u;)0L°V?i (s + - 13V2tan-\V2) - 8/n(6) + 16/n(z/; 
and so on. 

In the overlapping region, a(in)u{u) and k(^in)u{u) are obtained in the same way. As ^ 
approaches zero, a(^in)u{u) and k(^in)u{u) are expanded as 

ai^n)u{u) = -roCi'^ {{u -!)-(«- 1)^ + {u - if...) - vouCi'^ {{u -!)-(«- l)^...Xl26) 

- roiy'Cl'\u - 1)... + + u'A^J'> + z^^A^^^)... 

- V^r-ouHl''^ {{u -l)-{u- If...) - V3roiy'Hf\u - 1)... 



%„),(n) = iyKi'^ + u{u-l)Ki'^ + i^^Ki^^ + iy\u-l)Kf^ + iy'Kf^ (127) 

+ z/=^(w - l)Kf^ + v'K^^l.. - 2V?,rl{u - 1) {&f^ + z^C^^) + v^C^f^ + z^^Cf ...) 

+ 4r^z/i?« (3(ii - 1)2 + 2{u - If..) + 12ry (m - lfH^f\.. + i/^'i^'^^M - 1)... 

+ \rl 12 duu{u - .;)0L°)0i^ (^(« - l)^/n - 2(n - 1) 



- -{u-lf. 

In the overlapping region, the outer region solutions ay{u) and ku{u) in momentum space are 
given by 



= ^(^) y ^^^^(^-c^)CVr-^.(^-M«-l) + 3(^-1)^^-1) (128) 



«-l) +...) + O(ti-l)^ 



+ i (^-10 + 2y2tan"^(y2) -y27r + 4/n(6)' 
+ ^ ^40 - 1372^^^1(72) - 8/n(6) + 
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1 / / R 

k,{u) = -rl duu{p-u)(P^^^<p'^X(^{n-lfln{u-l)-l + V2tan-\V2) (129) 

. ,V2*».-HV2, - 3„. - 1) + (1^ - H _ 1^!^ 

- (" - 1)'] + '"tr + r»*r(« - 1) + o(u - if, 

where again is the Fourier transform of kiiy). We compare Eq fll26p . Eqf ll27p with Eq fll28p . 
Eq (fT^ respectively to decide that = = Cf'^ = A^J^ = Ri^^ = 0, 

+ 4/n(6) - 8ln{u)) 
and z/2i?(2) _ 2 / ^^^^^ _ ,^)0(o)0W^ (V2tan-\V2) - l) + rofc^ 

A^\k^^^ and Ri''^ are determined by matching with higher orders in the outer region solutions. 
With the coefficients determined in this fashion, two solutions are connected smoothly in the 
matching region. 

The constraint equations in the inner region are the same with those in the outer region. 
Plugging Ci^'^ = into Eq([Il9l), we obtain eI'^^ = 0. By using Cl^'^ = Ri^^ = 0, EqflT20D 
becomes 

u'Ei'^ = -—- / duuiu - (130) 



They are the same with the momentum space expression of Eq (!35|) and Eq (!36|l . 
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Appendix 

A Leading Corrections of the Toy-Model 
A.l Equations of motions 

With Eq( l29|) . the leading order Einstein equations for h{r,v), a{r,v) and k{r,v) become 
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Wr. = l-^[-l2r%{r,v)-2r{r^ + reo- pl)h'{r,v) (132) 

— 2rk'(r, v) + Ar'^pQa'{r, v) + 2k{r, v) + r'^k"{r, f )) 
1 ^2 fri-r^ 



W„„ = -^|M(_i2r''A(r,i.)-2r(r'' + rE„-pJ);i'(r,iO (133) 

— 2rk'(r, v) + Ar'^pQa'{r, v) + 2A;(r, f ) + r'^k"(r, f )) 

— ^(^-'^o)' + ^(2rr||i?(t;) - 2por',Civ)) = 

= -^{6r^h{r,v) + k{r,v) +r^Uo{r)h'{r,v) (134) 

— rA;'(r, f ) + 2r'^pQa'{r, v))) = 

,where Wu = W^x = Wyy and dots and primes indicate derivatives with respect to v and r 
respectively. The gauge field equations are 

Y"" = -^{poh'{r,v) + 2ra'{r,v) + r^a"{r,v)) =0, (135) 

2 

y = ^Civ) = 0. (136) 

The other components of the Einstein equations and gauge field equations are zero. These are 
the leading order equations in the naive derivative expansion. This means that t;-derivatives 
on h{r,v), a{r,v) and k{r,v) are ignored. 

Eq fll36p shows that there is no dynamics for the charge density. By the initial conditions 
mentioned in Secl2l C{v) = 0. A particular combination of Einstein equations, WrvUQ{r,v) + 
W^y = 0, gives 

E{v) = ^{d^<P^o){v)r (137) 

This equation indicates that E{v) ~ 0{e). This justifies that E{v) in the metric factor U{r,v) 
is suppressed by e to produce the second order terms in Eq( |T9i) . We solve Eqf ll3ip . Eqf ll35l) . 
Eq fll34p to get Eq (!3T!) . Eq (l33|) .Eq (!32!) respectively. The other Einstein equations are satisfied 
with the solutions. 

B Outer Solution in Extremal Limit 
B.l Dilaton solution 

We start with Eq fl20l) . In the case of extremality, Eq fl20|) can have a form of 
<P{u,v) = 0(o)(.) + |M!^(i_A,(.))-^^ (138) 

+ ^^^^^ ((Ai(t^) - 7)V2tan-\^^) + 2{A^{v) + 2)ln{u^ + 2u + 3)^ + A2{v). 
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As M — 7- oo, a boundary condition that we demand for the dialton is (f){u, f )|"~°° = 0(o)(v). This 
boundary condition yields A2{v) = — i-^ii'^) ~ 7)- With this, (j){r,v) has an asymptotic 
behavior of 

= 0(o)(.) + |^^(l-Ai(.))-5^ (139) 



^ dv(t)(Q){v 



36ro 
+ O(u-l), 



(Ai(t;) - 7)y2tan-^(V2) + 2(Ai(t;) + 2)/n(6) + ^(7 - Ai(^;))^ 



as n — )■ 1, near the black brane horizon. 



B.2 Metric and Gauge Field Solution 

As we discussed in Sec J3.3[ the regularity condition of the dilaton field forces Ai(z;) = 1. In 
this subsection, we follow this. For the extremal limit, we set eo = 2rQ and po = V^'^o- Eq fl3T]) 
becomes 

,(,,,„) = ,-..(.„) _ ±_(W? ( eju^- 1°-"- 15) _ i3y^,,„-.(l±i;) (140) 

^ ' ' ^ 864 \(^u-l){u^ + 2u + 3) ^2 

+ 16/n(w - 1) - 8ln{u^ + 2u + 3)) . 
As M — 1, this can be expanded as 

M^,-) = (,-9(^+27^"("-^) 

+ —(8 - l?>V2tan-^{V2) - 8/n(6)) + 0{u - 1) ) . 
216 / 

The near horizon expansions of Eq( |33ll and Eqf l32p are also given by 

« = -8 - + 1 /n n-1 - - + 72 - -7 tan-^ — =- 142 

864 ro \ M u u ^2 

+ 4(- + l)/Ti(M2 + 2ti + 3) +a2(i;) v o iv ; 



A:K^;) = -^r2(9,0(o)(t;))'(^^(n -1)2(^2 + 2m + 3)/n(^i-l)-i/2 (143) 
_ y2tan-i(^) (^-« + ^(«- 1)2(^2 + 2^ + 3) 

- - 2n + 3) ^^^^, + 2. + 3) + 1(.2 - 10. - 15)(. - 1) 

27 18 

+ rQu\[v) - 2V3r2ai(t;) + 2r^(w - 1)2(^2 + 2m + 3)/ii(t;), 

The integration constants,/ii(f ), ai(f), a2{v) and A;i(t;) are determined by the boundary con- 
dition flM|l . As li — )■ 00, the asymptotic expansion of k{u,v) can have terms of O(u^). These 
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are non-normalizable modes which give deformation of the boundary metric. The terms are 
removed by imposing hi{v) = — y|^^ (5^,0(0) (t'))^. Near AdS^ boundary, a{u,v) presents 0(1) 

and 0{^) terms. The former corrects the chemical potential and the later does the charge 
density. To eliminate these terms, ai{v) and a2{v) should be properly chosen as di{v) = and 

Near the black brane horizon, the behavior of the leading back reactions are given by 
„.,„.„, . ^ (144) 

+ 7^(8 + - lsV2tan~\V2) - 8/n(6)) + 0{u - 1) 

21b 2 

a{u,v) = ^iM5lMl!(_8/n(ii-l)(3-2(«-l) + O(«-l)2)- 30 (145) 
864 To 

+ 6V2tan-\V2) - Sv^tt + 12/n(6) 

+ 1^40 - 13V2tan-'\V2) - 8/n(6) + j (n - 1) + 0{u - if 

= -\rl{d,ct)i^^){v)f {^{u - lfln{u - 1) - 1 + V2tan-\V2) (146) 

^ (yl,„„-.,v^)-3)(.-l,.(l^-H_f!^ 



/c(m, v) 



9^2 



{u - 1)M + roA;i(t;) + r^kiivfu - 1) + 0{u - 1)' 



C Equations in Extremal Backgrounds with z;-derivative 
Retained 

In this section, we develop the Einstein equations([5]) and the gauge field equations([6]) without 
ignoring f-derivatives. We start with Eq( l29|) . The only assumption in this section is that the 
equations are linear in h{u,v), a{u,v) and k{u,v). The Einstein equations are 

9h'(ii u) 1 

rlWrr = ^-^--ducPiu,v)du<Piu,v) = 0, (147) 

u 2 

1 / 2u'^ ■ 

= —r{-12u%(u,v) + 2u(3-u^-2u)h'(u,v) + h'(u,v) (148) 

2u^ \ To 

+ a'{u, f ) + 4: (2A;(m, v) - 2uk'{u, v) + u^k"{u, vf) , 

ro ' 

- 7^dy(l){u, v)du(j){u, v) = 0, 
2ro 
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W ^ W..+ { ' ) (149) 



""^ ' 2VlC{v) + 2uE{v) - 2{u^ -Au + 3)h{u, v) + 



^dy(l){u, v)d^(j){u, v) - ro (- — ^ ) (9„(/)(m, v)dv4>{u, v) = 0, 



^2 



-^Wii = —6u'^h(u,v) — u(u^ — 4:U + 3)h'(u,v) — a(u,v) (150) 

ro 

+ ^ (-uA;'(M,f) - A;('u,i;)) = 0, 

where u = the rescaled radial coordinate again. The primes and dots denote derivatives 
with respect to u and v respectively. Gauge field equations are 

2 2 

L^Y' = Civ) + V3hiu, v) + —a'{u, v) = 0, (151) 
To ro 

r'^Y" = -2ua'{u, v) - V3roh'{u, v) - u'^a"{u, v) = 0. (152) 

As discussed in the beginning of SecEl for introducing the scaling( l42i) the equations in 
position space of v should be transformed to the momentum space. It is worth showing how 
one gets an expression Eq (ll47l) in momentum space, for example. We define the Fourier 
transform of h{u,v) by 

/oo 
hMe^'^^'du. (153) 
'OO 

Substituting of Eq(|59l) and Eqf ll53p into Eq fll47l) and acting an integral operator ^ e'^'^^dv 
on it, we get 

2h' (u) 1 f°° 

rlw„ = ^ - - / duduCpu^MduMu) = 0. (154) 

To deal with other equations, we define Fourier transforms of a{u, v) and k{u, v) as that of 
h{u,v). For E{v) and C{v), 



E{v) = / e'^'^'E^du, (155) 

/oo 
-oo 

With these, the other Einstein equations in momentum space are given by 



1 f 2u'^ 

— -Uu^hju) + 2u(3 -u^- 2u)h'iu) + iv h'iu) (156) 

2u^ \ ro 



+ ^^^<(^/) + ^ [2K{u) - 2ukl{u) + u'^k'liu)) 
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— / duji{u - Lo)(j)^^^{u)du(f)uj{u) = 0, 
Zr 



4 _ 4rr^ + 3rn 



r 



W = W,,+ [ ] Wr. (157) 



- 4m + 3 
2m 



~ ^0 ( —2 ) / dui{p ~ uj)(j)^^^{u)du<i)u^{u) 



oo 



+ ^1 du{u - u)uJ(l)u-oj{u)(t)^{u) = 0, 

^ J -co 

r2 2^/3?/2 

'-Wii = -Qu^K(u)-u(u^-Au + 3)h'(u)-=-^a'(u) (15^ 
ro To 

+ — {uk'^{u) - K{u)) = 0. 



The gauge field equations become 



.2 



= (c, + VSKiu) + ^<(m)^ = 0, (159) 



^2yv ^ -2ual{u) - v^rohliu) - m^<(m) = 0. (160) 

D Counting Power of e 

In this section, we argue the parametric order of the inner region solutions in e. The basic idea 
is that we transfrom the inner region solutions to the position space and check their powers of 
e. For a simple example, we discuss the dilaton solution. We design the zeroth order dilaton 
solution in position space as in Eq ([TT|) . By Fourier transformation, we obtain its expression in 
frequency space as 

^(0) -J_r e-'^-f{—)dv. (161) 



Scaling of the integration variable v a.s t = ev takes this expression to 

~ r e-^Trfi^L)dr^-g(^\ (162) 
2718 J _^ ro e e 

where g{-) becomes an 0(1) function. By observing Eqf ll62l) and dilaton solution in the inner 
region fBSj) . one can recognize that each subleading correction to the dialton field in momentum 
space can be written as 

= \9{^) {h(r)iO + , (163) 

where (p^f^^^^^ denotes ith order correction in small frequency to the dilaton field. /i(j)(0 and 
a(j)(^) are functions of C, only. It turns out that the terms multiplying a{i){C,) produce terms 
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which are proportional to e^ln{e) in the position space. It is obscure to count power of e of 
these terms. In this discussion, we exclude these. The perturbation expansion of the dilaton 
solution becomes a form of 

(pi^n). = ^^(^) (1 + + ^'h2){0-) , (164) 

up to the logarithmic terms. Fourier transformation defined in Eqf l59l) takes this expression to 
position space, which is given by 

e^^^du-g{-) (1 + a;/i(i)(0 + ^'/i(2)(0-) ■ (165) 

-oo £ S 

Again we rescale the integration variable u a.s u = eu, then the expression becomes 

/oo 
e^'^'^'dugiu) (1 + + e^u^h^2)iO-) (166) 

-oo 

where F^'^^ev) = e^^^'"dcdg{uj) and the prime indicates derivative with respect to its argu- 
ment. The property of F^^\ev) as noted in Eq f lT^ shows that F^^\ev) and its derivatives 
with its argument are 0(1) functions. Compare Eq( ll64p to Eqf ll66p . This shows that counting 
power of u in the momentum space is the same as the counting power of e in the position 
space. 

We apply this argument to the back reactions in the inner region. For the simplest case, 
let us check h^(^n^^{^) with Eqf l95p . This contains derivative of the dilaton field, which can 

be expressed as <P^(^n)uj{j;0 = e/(f)^(^0- We switch Eq fl95l) to position space with Fourier 
transformation as defined in Eqf ll53p and rescale integration variables as in the discussion of 
the dilaton. This time, we scale u as well as u in Eq fl95p . Then, we obtain following expression: 

h%{v,0 = e r e^'^^dudcoco\u-oormf{u-u)g{^09i'^0 (167) 

J J —oo 

= eG{ev), 

where G{ev) is an 0(1) function. Consequently, it turns out that h^(2i)i'^-:0 is in the first order 
in e in position space. Comparing Eq fll67p to Eq fl^?^ . one can recognize that the correct power 
counting of the small frequency in the momentum space is to count not only u but also oj and 
— a; in the integrand of Eq (!95|) . By the the same argument, we get 

~ C(^'0~O(£), (168) 

- «M(^'0~fcJ2)(-,0-O(£^), 
«S)(^,0 - k\^Uv,0-O{e') and fcgj., ~ 0(.^). 



30 



References 

[I] Sayantani Bhattacharyya, Veronika E Hubeny, Shiraz Minwalla, Mukund Rangamani, 
JHEP 0802:045(2008), [arXiv:hep-th/0712.2456] 

[2] Sayantani Bhattacharyya, R. Loganayagam, Shiraz Minwalla, Suresh Nampuri, Sandip P. 
Trivedi, Spenta R. Wadia, JHEP 0902:018 (2009), [arXiv:hep-th/0806.0006] 

[3] James Hansen, Per Kraus JHEP 0904:048 (2009), [arXiv:hep-th/081 1.3468] 

[4] J. Erdmenger, M. Haack, M. Kaminski, A. Yarom, JHEP 0901:055 (2009), [arXiv:hep- 
th/0809.2488]; N. Banerjee, J. Bhattacharya, S. Bhattacharyya, S. Dutta, R. Loganayagam, 
P. Surowka, [arXiv:hep-th/0809.2596]; 

[5] Giuseppe P ohcastro ,Dam T. Son, Andrei O. Starinets, JHEP 0209:043(2002), 
|hep-th/0205052|, Giuseppe Po licastro ,Dam T. Son, Andrei O. Starinets, JHEP 
0212:054(2002) |hep-th /0210220| 

[6] R. Baier, P. Romatschke, Dam T. Son, Andrei O. Starinets, M. A. Stephanov, JHEP 
0804:100(2008) [arXiv:hep-th/0712.2451] 

[7] Sayantani Bhattacharyya, Subhaneil Lahiri, R Loganayagam, Shiraz Minwalla, JHEP 
0809:054(2008) , [arXiv:hep-th/0708. 1770] 

[8] Dam T. Son, Andrei O. Starinets, Ann. Rev. Nucl. Part. Sci. 57:95-118(2007), [arXiv:hep- 
th/0704.0240] 

[9] Giuseppe Policastro ,Dam T. Son, Andrei O. Starinets, Phys. Rev. Lett. 87, 081601 (2001), 
[arXiv:hep-th/0104066j 

[10] Makoto Natsuume, Takashi Okamura, Phys. Rev. D77, 066014(2008), [arXiv:hep- 
th/0712.2916] 

[II] Thomas Faulkner, Hong Liu, John McGreevy, David Vegh, [arXiv:hep-th/0907.2694] 

[12] Mohammad Edalati, Juan I. Jottar, Robert G. Leigh, JHEP 1001:018(2010) [arXiv:hep- 
th/0910.0645] 

[13] Mohammad Edalati, Juan I. Jottar, Robert G. Leigh, JHEP 1004:075(2010) [arXiv:hep- 
th/1001.0779] 

[14] Adel Awad, Sumit R. Das, Archisman Ghosh, Jae-Hyuk Oh, Sandip Trivedi, Phys. Rev. 
D80, 126011 (2009). [ar Xiv:0906.3275 [hep-th]]. 

[15] M. J. Duff, J. T. Liu, Nucl. Phys. B554(1999) 237-253, (hep-th/9901149l; D. Klemm, 
Nucl.Phys. B545 (1999) 461-478, [ hep-th/9810090j 



31 



